Point symmetries in the Hartree-Fock approach: Densities, shapes and currents 
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Three mutually perpendicular symmetry axes of the second order, inversion, and time reversal can 
be used to construct a double point group denoted by D™. Properties of this group are analyzed 
in relation to the symmetry and symmetry-breaking effects within the mean-field (Hartree-Fock) 
theories, both in even and odd fermion systems. We enumerate space symmetries of local one-body 
densities, and symmetries of electromagnetic moments, that appear when some or all of the D^iP 
elements represent self-consistent mean-field symmetries. 
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I. INTRODUCTION 

Point-group symmetries play a very important role 
in nuclear mean-field theories. Two distinct aspects of 
this role can be singled out. Firstly, the point symme- 
tries of a Hamiltonian provide good quantum numbers 
that can conveniently be used to label its eigenstates. 
They help to formulate selection rules for electromag- 
netic and/or other types of transitions, and allow for 
solving the stationary problems in subspaces rather than 
in the complete Hilbert space of the problem in question. 
In that respect the use of point symmetries in nuclear 
physics resembles their use in other branches of physics. 
Secondly, however, and this aspect is more specific to 
nuclear structure domain, the use of the self-consistent 
Hartree-Fock (HF) or Hartree-Fock-Bogolyubov (HFB) 
mean field approximations invariably leads to the prob- 
lem of self-consistent symmetries and related sponta- 
neous symmetry-breaking mechanisms |l| . 

In this article we aim at describing properties of nu- 
clear one-body densities under the action of point sym- 
metries. For the time-even densities we calculate the 
electric multipole moments which give information about 
nuclear shapes. Various point symmetries obeyed by the 
hamiltonian lead then to various types of allowed shapes. 
In addition, for the time-odd densities we calculate mag- 
netic multipole moments which give information about 
current distributions in nuclei, i.e., about the "shapes" 
of matter flow. Again, various conserved symmetries re- 
strict these flow patterns in different ways that are stud- 
ied in this paper. 

Numerous experiments indicate that a great number 
of nuclei are deformed in their ground states. Interpreta- 
tion of the corresponding results shows that most often 
the shapes involved are axially-symmctric. Many realis- 
tic calculations, e.g., those based on the nuclear mean- 
field approximation, reproduce these experimental data. 
However, the same calculations suggest that the excited 
nuclear states often correspond to the nucleonic mass dis- 
tributions that have the so-called triaxial shapes. It thus 
becomes clear that in a realistic description of the nuclear 
properties, the spontaneous symmetry breaking leading 



to the triaxially symmetric objects must be given atten- 
tion. 

The classical point group that contains three mutu- 
ally perpendicular symmetry axes of the second order 
passing through a common point is denoted D2 (cf., e.g., 
Refs. [@-|5|]). Attaching to D 2 the three mutually perpen- 
dicular symmetry planes spanned on the symmetry axes 
gives the D 2 h point group which contains all the spatial 
symmetries of interest in the present paper. 

As it is well known, the classical {single) point groups 
can be applied to spinless particles and/or systems of 
an even number of fermions, and thus to even-even nu- 
clei. However, for odd fermion systems and, in partic- 
ular, in the single-nucleon space, these have no faithful 
irreducible representations. There exist two methods to 
remedy this problem. One is to extend the notion of the 
group representation and to introduce projective or ray 
representations ■ Another one, which is employed 

in the present work, is to enlarge the single groups by 
adjoining the rotation through angle 2tt and all its prod- 
ucts with the original group elements, and to double in 
this way the order of the group |^] . 

Physically, the need of such an extension is related to 
the fact that in the space of spinors the rotation through 
angle of 2ir necessarily changes the sign of the wave func- 
tion of an odd-fermion system. Since within the group 
theory a multiplication of a group element by a number 
is not defined, the change of sign must be introduced as 
an extra group element. The point group enlarged in this 
way is called the double point group and usually denoted 
with the superscript "D" (cf., e.g., Refs. UH), although 
some authors, see, e.g., [||, denote single and double point 
groups by the same symbols. Here we follow the former 
convention, and thus the double group corresponding to 
single group D 2 h is denoted by D? h . 

In the case of classical objects the elements of a symme- 
try point group are real orthogonal coordinate transfor- 
mations. In quantum mechanics it is often of advantage 
to take into consideration both the spatial symmetries 
and the time-reversal operator explicitly and treat them 
as elements of a common ensemble of symmetry opera- 
tors. Time-reversal symmetry operator (antilinear) and 
the space symmetries (linear) have usually been consid- 



1 



ered separately (cf., e.g., Ref. ||). Here we follow Ref. 
H and add the time-reversal operator to the set of group 
elements, thus obtaining new groups denoted Dj h and 
D^. Hence, the D 2 h group is a subgroup of Dj h com- 
posed of its linear elements, and similarly, the DJ? h group 
is composed of the linear elements of D ™ . 

Gauge symmetries, which pertain to pairing correla- 
tions of nucleons, can be added independently and are not 
considered in the present study. In particular, neutron- 
proton correlations are not discussed. Therefore, the 
isospin degree of freedom is irrelevant in the discussion 
and can be disregarded to simplify the notation. 

In this paper our goal is threefold: First, in Sec. O, 
we present and discuss properties of the single group D^ h 
and double group D™ that are appropriate for a descrip- 
tion of even and odd fermion systems, respectively. In 
particular, we recall the classification of representations 
of both groups, and classify properties of the group ele- 
ments when the y a re represented in fermion Fock space. 
Second, in Sec. 
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we present explicit symmetry prop- 
erties of local densities with respect to the operators in 
the D™ group. This problem has been solved in particu- 
lar applications PlJlOfl; however, it can be solved in many 
different ways, and it is useful to have a systematic ap- 
proach which enumerates all available options. Although 
the local densities are most important for applications 
using the local density approximation (LDA), or those 
using the Skyrme effective interaction (see respectively 
Refs. plH or Ref. Q for reviews), they also define 
general properties of avera ge v alues of any local one-body 
operator. Finally, in Sec. [IVL we discuss symmetries of 
multipole moments which define the nuclear shapes and 
currents, and in Sec. ^ we present conclusions which can 
be drawn from our study. In the companion paper [pr[ , 
we discuss physical aspects of the symmetry-breaking 
schemes pertaining to the point groups in question. 



II. SYMMETRY OPERATORS 

The point groups of interest in this paper can be in- 
troduced in two ways. The first one consists in defining 
an abstract point group by giving its table of multipli- 
cations, and then classifying the states and operators in 
the fermion Fock space according to the relevant irre- 
ducible representations (irreps). The second one, which 
we follow below, is more intuitive, and consist in defin- 
ing the symmetry operators in the Fock space first, and 
then identifying their multiplication tables and the cor- 
responding group structures. 



A. Fock-space representations 

It will be convenient to use the Cartesian representa- 
tion of the symmetry operators. Let I k for k = x,y, z 
denote the Cartesian components of the total angular 



momentum operator (generators of the group of rota- 
tions). In the coordinate-space representation these op- 
erators read (^=1) 



n—l n—1 



iV + 



2"k 



(1) 



where jj^, 



and |cr fe 



(«) 



are operators of the total, 



orbital, and intrinsic angular momenta, respectively, of 
particle number n. By definition, operators (|lj) act in the 
Hilbert space Ha of A-particle states, and the number of 
particles A appears explicitly in their definitions. 

One can use another representations of I k , the so-called 
second-quantized, or Fock-space form, 



(2) 



where (At|jfc|^) are the matrix elements of the angular- 
momentum operators in the single-particle basis defined 
by the fermion creation and annihilation operators a+ 

and a u . Operators I k in the form of Eq. (|2|) do not ex- 
plicitly depend on A, and act simultaneously in all the 
A-particle spaces, i.e., they act in the Fock space H, 



H = Hq © Hi 



(3) 



In each subspace Ha, operators (g) are equal to (Q). 
Since both act in different domains, one should, in prin- 
ciple, denote them with different symbols. However, one 
usually understands definition ([!]) as a prescription to 
construct I k for all values of A simultaneously (adjoined 
by lk=0 for A—0). With this extension, operators (|l|) 
and (|^) are equal. In this section we understand that 
all operators act in the Fock space (||), while the cor- 
responding definitions are given in the coordinate-space 
representation. 

We introduce three standard transformations of rota- 
tion around three mutually perpendicular axes, Ox, Oy 
and Oz, through the angles of ir as 



R k ee e-™ 1 " = 



(4) 



Similarly, we introduce three operators of reflection in 
planes yz, zx, and xy, for k — x,y, z, respectively, which 
can be written as 



S k = PR k , 



(5) 



where the inversion operator is denoted by P. The order 
of operators in Eq. (o) is unimportant because 



[P,R k ] = 0. 



(0) 



Finally, the (antilinear) time-reversal operator in the 
coordinate-space representation is defined as JTa : 



2 



K. 



(7) 



where K is the complex conjugation operator associated 
with the coordinate representation. 

In what follows it will be convenient to denote with 
separate symbols the products of T with P, Pfc, and Sk 
iflql , i.e., the seven additional (apart from T itself) anti- 
linear operators read 



pT 



PT, 

PfcP, 

S k f. 



(8a) 
(8b) 
(8c) 



The order of multiplications in the above definitions is 
irrelevant since 



[P ) f] = [R k ,f} = [S k ,f} = 0. 



(9) 



In nuclear physics applications the linear operators P, 
Rk, and Sk are usually referred to as inversion, signature, 
and simplex. The antilincar operators P T , R^, and S J 
will be called T-inversion, T-signature, and T-simplex, 
respectively. 

For completeness, yet two another operators must be 
added to the above symmetry operators. One of them is, 
of course, the identity operator, E, which can be treated 
as the rotation through angle equal to or 4tt about an 
arbitrary axis. The second one is the rotation through 
angle 2n about an arbitrary axis, i.e., 



E 



-i2nLk 



(-^)=(-l) A E (10) 



where &o is the unity 2x2 matrix. We see that only 
for even systems E is equal to identity, while for odd 
systems it is equal to the minus identity. We should 
keep in mind, that in the group theory there is no such 
a notion as a change of sign. Operators like (— 1) A may 
appear in representations, like here they do appear in the 
Fock-space representation, however, one can not u se t hem 
when defining the group structures in Sees. [IB and [IC 
below. 

To investigate multiplication rules of the symmetry op- 
erators introduced above one explicitly calculates prod- 
ucts of them. For example, the products of two signatures 
are 



n=l 



/ - (n) - 



(n) 



(11) 



and the square of the time reversal reads 

A 



n=l 



(-(*^f)=E. 



(12) 



It is obvious that these results depend on whether A is 
even or odd. 

Therefore, in what follows we introduce notation which 
explicates whether the operators act in even or odd 
fermion spaces, 7i+ or 7i- , 



Ho 
Hi 



• n 2 
•n 3 



T~(-A=2p+1 



(13a) 
(13b) 



Any Fock-space operator U:H — >H, which conserves the 
particle number, is split into two parts with different do- 
mains, i.e., 



U = U + U, 



(14) 



where the bold symbols denote operators which act in 
the even- A spaces, while the script symbols denote those 
acting in the odd- A spaces, i.e., 



U:H+ — ►«+, 
U:H- — > H- 



(15a) 
(15b) 



With these definitions we are now in a position to inves- 
tigate the group structures appearing for the introduced 
operators. 



B. Single group D 2h for even systems 

For even fermion numb ers w e consider the Fock-space 
operators defined in Sec. [I A, and restrict them to 
Eqs. ( |T4] ) and (15a). Then, the complete multiplication 
table reads 





Kl = 




= E, 


(16a) 






) =P 2 


= E, 


(16b) 




RfeSfe 


= SfcRfc 


= P, 


(16c) 






= s fe Rfe 


= P, 


(16d) 


RfcR-fc — Rfe R* 


= s fc s fc 


= s fe s fc 


= t, 


(16e) 


Rfe Sfe = Sfe Rfe : 


= s fc 


= SfeRfe 


= pT 


(16f) 


= x, y, z, and 










RfeR; = SfcS/ = 


R fe Rf : 


= slsf -. 


— Rm j 


(17a) 


RfeS; = SfeR; = 


Rfc s f = 


z Sfe : 


— S m , 


(17b) 


RfeRf = Rfc R/ 




= S^S; = 


- R T 

""m > 


(17c) 


Rfe S; = SfeR^ = 


= RfcSf 


= S^R; = 


_ cT 

~ °m> 


(17d) 



for (k,l,m) being an arbitrary permutation of (x, y,z). 

We see that the 16 operators acting in the even-A 
fermion spaces constitute the Abelian single group which 
we denote by Dj h , 



{E, P, R fe , Sfe, T, P T , Rfe , Sfe }, 



(18) 
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for k — x,y, z. The half of the elements in Eq. ( |l8| ) are 
linear operators and the other half are antilinear. 

It follows from the multiplication table of the Dj h op- 
erators, Eqs. ( |l6|) and (|l7f), that the whole group can 
be generated by its four elements only. These elements 
are called the group generators. Various possibilities of 
choosing the generators are discussed in Ref. here 
we only mention that, e.g., the subset {T, P, R x , R y } 
can be used to obtain all the operators that belong to 

The Dj h group has two important subgroups, the 
eight-element Abelian group D 2 h composed of all the lin- 
ear Dj h operators, 



D 2h : {E,P,R fc ,S fe }, 
and the four-element Abelian group D2, 
D 2 : {E,R fc }. 



(19) 



(20) 



Obviously, the D 2 h subgroup of D J h , being an Abelian 
group of order eight, has eight one-dimensional irreps. 
These can be labelled by eigenvalues equal to either +1 
or —1 of three of its generators, say, P, R x and Ry. 

We introduce names for the (one-dimensional) bases 
of the associated irreps according to the following con- 
vention. First, a basis is called invariant if it remains 
unchanged (belongs to eigenvalue +1) under all three sig- 
nature operators R&, and it is called either x-, or y-, or 
z-covariant if it transforms under the signature operators 
like the x, or y, or z coordinates, respectively. Secondly, 
prefix pseudo is added for bases which are odd, i.e., be- 
long to eigenvalue —1 with respect to the inversion P. 

Since T is an antilinear operator and also an invo lutiv e 
operator [i.e., its square is equal to identity, Eq. (16a)], 
we can always choose the phases of all the basis states 
so that t hey b elong to the eigenvalue T =+1 of T JlH] 
(see Sec. [IE below). In this way we construct eight 
irreducible one-dimensional corepresentations (ircoreps) 
of Dj h , all being even with respect to the time reversal. 
After Wigner |1, we call the representations of a group 
containing antilinear operators corepresentations to em- 
phasize the fact that they are not the representations in 
the usual sense (see Appendix for details). By a suitable 
change of phases of the basis states we can obtain an- 
other set of eight ircoreps of Dj h , all of them odd (i.e., 
belonging to the eigenvalue T=— 1) with respect to the 
time reversal. We use prefix anti to name these time-odd 
ircoreps. Obviously, time-even and time-odd ircoreps are 
pairwise equivalent (cf. Ref. ||). 

Note that all operators acting in the even fermion 
spaces H+ can also be classified according to the same set 
of sixteen ircoreps of Dj h . All these ircoreps are listed in 



Table [II together with explicit transformation properties 
of several examples of one-particle operators belonging to 
each ircorep. 



C. Double group DjjP for odd systems 

For odd fermion numb ers w e consider the Fock-space 
operators define d in Sec. [I A , and restrict them to 7i_, 
Eqs. @ and fll5b| ). Since operator £ [odd-fermion- 
number part of E of Eq. (|l0|)] is now an independent 
group element, additional partner operators should be 
introduced in order to construct the double group D™, 
i.e, 



P = £V, r = £T, V 1 



EV 1 



(21) 



and 



K k = £K k , S k = £S k , Til = Si = ESl, (22) 

for k=x,y,z. Now the group multiplication table reads 



for k 









si = f 2 


= £, 


(23a) 








) =v 2 


= £, 


(23b) 






n k s k 


— SkTZk 


= v, 


(23c) 








= sln T k 


= v, 


(23d) 






= SkSl 


= SlSk 


= f, 


(23e) 


K k Sl 


= sln k 


= TZjSk 


= Skl^l 


= P T , 


(23f) 


= x, y, z, and 










n k iii 


- S k Si — 


K T k nJ-- 


-slsf- 


- V 


(24a) 


■RkSi = 


-- s k iii = 


ais? = 


-- slnj -- 


— ^rn 1 


(24b) 


n k fif 




= SkSf 


= SlS[ -- 


= a T 

lx "m 1 


(24c) 




= S k TZf -- 


-- HkSj 


= sl^ -- 




(24d) 



for (k,l,m) being an odd permutation of [x, y, z), while 
relations identical to ( ]l7| ) hold for an even permutation 
of (x, y, z). After multiplying relations fl23| ) and ( p4] ) by 
£ once or twice, one can easily obtain the remaining el- 
ements of the multiplication table, i.e., those which per- 
tain to products involving one or two partner operators 



The D™ group is thus composed of 32 operators: 
D™ : {£,V,T,f> T ,K k J k ,KlJl, (25) 

£,P,T,P T ,1i k ,s k ,1il,sl}. 

One can see that this double group is not Abelian, be- 
cause relations (|24|) now do depend on whether the per- 
mutation (k, I, m) of (x, y, z) is even or odd, whereas for 
the single group, relations (|17|) are independent of that. 

On e m ay note that the Fock-space operators U of 
Sec. |ll A| and the odd-fermion-space operators U of 
Eqs. (|14|) and ( |l5b| ) obey exactly the same multiplication 
rules of the double D™ group. Therefore, one might, 
in principle, consider only the double group D™ and 
refrain from studying the group structures in even and 
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odd spaces separately. However, at the level of repre- 
sentations, one would then have been deprived of impor- 
tant properties of operators like T 2 =E or T 2 =—£ (see 
Sec. [IE ), neither of which holds in the whole Fock space, 
cf. Eqs. |L2|) and @. 

Since the squ ares of the time reversal, signatures, and 
simplexes, Eq. (23a), are equal to E, the whole double 
group D™ can be generated by the same operators which 
generate the single group in even systems, Sec. II B. So 
the double group also needs four generators; for instance, 
the set of four elements, T, P, 1Z X , and IZy, can be used 
to obtain the entire double group of the DjjP operators. 

The linear operators of the D™ double group form the 
sixteen-element double group D; 



2h> 



u 2h 



{£,V,K k ,S k ,£,P,K k ,S k }. 



(26) 



This group has 10 equivalence classes (cf. Refs. [2,3,5]). 
There are 6 classes composed of 2 elements each, i.e., 
{P k ,P-k} and {5^,^} for k = x,y,z, while the remain- 
ing elements: {£}, {£}, {P} and {P} form 4 one-element 
classes by themselves. The group is not Abelian and pos- 
sesses, apart from the 8 one-dimensional irreps already 
known for the single group, another 2 two-dimensional 
spinor irreps. The spinor irreps can be labelled by the 
parities (the eigenvalues 7r=+l or —1 of the inversion 
operators V), see Appendix. 

The time reversal T is an antilinear and also an antiin- 
volutive operator [i.e., its representations give the minus 
identity when squared, Eq. fl23a|)1 , and therefore it cannot 
be diagonalized |l5| (see Sec. [II E| below), and used for la- 
belling the ircoreps. A Hermitian antilinear involu- 
tive oper ator, i.e., either a T— signature or a T— simplex 
[see Eq. ( [23b| )l should be chosen to serve this purpose. 
For instance, a pair of commuting Hermitian operators 
P and V,y can be used to label the D™ ircoreps. As for 
the single group Dj h , ircoreps being either even or odd 
with respect to P y can be obtained one from another 
by a suitable change of phase, and are therefore equiv- 
alent. In analogy to the one-dimensional ircoreps, the 
bases of spinor ircoreps belonging to pairs of eigenvalues 
of {P,P-T} equal to {+1,+1}, {-1,+1}, {+1,-1} and 
{ — 1, —1} can be called the spinor, pseudospinor, anti- 
spinor and antipseudospinor bases, respectively. 

Note that only spinor ircoreps of D™ appear in the 
classification of states of systems with odd numbers of 
fermions. However, the operators acting in 7i_ can all 
be classified according to the one-dimensional ircoreps 
of D™, similarly as operators acting in 7i + can all be 
classified according to the corresponding one-dimensional 
ircoreps of Dj h . (This is completely analogous to the 
fact that fermion-number conserving operators can carry 
only integer angular momenta, i.e., they are integer-rank 
tensors.) Therefore, whenever we consider the action of 
the Dj h or D™ operators on fermion states we always 
specify whether they act in even TL + or odd W_ spaces, 
and use for them the corresponding notations U and IA 



of Eqs. (15a) and (15b). On the other hand, whenever 
we consider transformation properties U^OU of operators 
O with respect to the D™ group, we do not make this 



distinction, and use for them notation U of Sec. II A 



D. Cartesian harmonic oscillator basis 

One often uses the Cartesian harmonic oscillator (HO) 
basis to solve the self-consistent equations when nei- 
ther spherical nor axial symmetry is assumed, see, e.g., 
Refs. |^|l(J. The Cartesian HO states are identified 
by the numbers of oscillator quanta, n Xl n y , and n z , 
in the three Cartesian directions, and by the spin pro- 
jection s z =±i on the z axis. For the standard HO 

phase convention, this basis is real, K,\n x n y n z , s z ) = 
s 2 ), where according to our standard conven- 
tion the script symbol K, denotes the coordinate-space 
complex-conjugation operator acting in the odd-fcrmion- 



number space 7i.~, see Sec. II A. 
For the HO states the fol 



owing relations hold 



V\n x n y n z , s z ) = (-1) 
f\n x n y n z , s z ) = (-1)^" 
K x \n x n y n z , s z ) = i(-l)™ 5 ' 
P y \n x n y n z ,s z )= (-1) 



)!,„+ri„+n z 



z \n x n y n z ,s z ), (27a) 
i x n y n Zl -s z ), (27b) 
¥1 \n x n y n z , -s z ), (27c) 
" - +i - s *\n x n y n z ,~s z ), (27d) 
\n x n y n z ,s z ), (27e) 
(27f) 



K z \n x n y n z , s z ) = 
£\n x n y n z ,s z ) = -\n x n y n Zl s z ), 

from where one can find similar equations for all the re- 
maining operators of group D™. Since the HO Hamil- 
tonian is symmetric under D™, its eigenstates can be 
classified according to the ircoreps of D™. It is easily 
seen that the HO states \n x n y n z , s z ) form bases of the 
spinor, pseudospinor, antispinor and antipseudospinor 



ircoreps for {N=n x 



-n z , N y =n x +n z } being {even, 



odd}, {odd, odd}, {even, even} and {odd, even}, respec- 
tively (see Appendix). The entire HO basis would have 
belonged to the spinor and pseudospinor ircoreps only, if 
the basis states and phase convention were chosen differ- 
ently, see Ref . [H . 



E. Properties of the Dj h and D™ operators 



In this section we recall properties of the Dj h and 



T)TD 



operators when they are represented in the fermion 
Fock space. Within representations, apart from the 
corresponding multiplication tables, Eqs. (p~6j) — (JT7|) and 
(p3|)-(p^), these operators are characterized by their 
hermitian-conjugation properties. Since all the Fock- 
space representations of the Dj h and D™ operators are 
unitary, they are hermitian or antihermitian depending 
on whether they are involutive or antiinvolutive, respec- 
tively. Properties of these operators are very different de- 
pending on whether they arc linear or antilinear. These 
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characteristics are summarized in Tables § and O, where 
the and D^j operators are split into two or four 
subsets, respectively. Below we review the properties of 
operators in each such subset. 

TABLE I. Properties of the Djh operators U in even 
fermion spaces. 





Linear 


Antilinear 


Hermitian (U 2 =E) 


P, R-fe, Sfc 


RJ, SI, T, P T 



TABLE II. Properties of the D™ operators U in odd 
fermion spaces. 





Linear 


Antilinear 


Hermitian (U 2 —£) 


V 




Antihermitian (U 2 ——£) 


TZk, S k 


f, V T 



For each linear Dj h or D™ operator one can attribute 
quantum numbers to fermion states. These quantum 
numbers can be equal to ±1 or ±i for hermitian (in- 
volutive) or antihermitian (antiinvolutive) operators, re- 
spectively. Therefore, the parity operators P or V give 
the parity quantum numbers, 7T=±1, the signature oper- 
ators Rfe give the signature quantum numbers, r=±l, in 
even systems and the signature operators TZk give r=±i 
in odd systems. Likewise, the simplex operators Sfe and 
Sk give the simplex quantum numbers, s=±l and s=±i, 
respectively. 

Antilinear operators do not give good quantum num- 
bers, and their role is very different, depending on 
whether they are hermitian or antihermitian, Tables || 
and[T§ 

For each hermitian antilinear Dj h or D™ operator, 
i.e., for R£, Sj, T, P t , TZ\, or one can find a 
basis consisting solely of its eigenstates with the com- 
mon eigenvalue equal to 1 |l5| . Indeed, if state \9) is 
an eigenstate of, e.g., Rj, the corresponding eigenvalue 
must be a phase, i.e., R^\^)=e 2t< l'\^) . In such a case, 
state |\E' / )=e l ^| v E'} is an eigenstate of R ]» with eigenvalue 
1. This demonstrates explicitly that properties of eigen- 
states of are, of course, phase-dependent. In the case 
when state |^) is not an eigenstate of Rj, one can trans- 
form the two linearly-independent states \$f) and R 



into eigenstates of it£ with eigenvalue 1 by symmetriza- 
tion and antisymmctrization of the two: 

|*.)= |*)+ R%\*), 
\* a ) = i\*)-iRZ\*), 



(28a) 
(28b) 

which also requires a specific phase. Therefore, phase- 
convention properties of states are essential for a discus- 



sion of bases of eigenstates of the hermitian antilinear 
operators, and in Ref. |Q a special discussion is devoted 
to this problem. 

One should also remember, that only linear combina- 
tions of basis states with real coefficients remain eigen- 
states of any hermitian antilinear operator. This is in 
contrast to properties of linear operators, for which a 
linear combination of eigenstates, corresponding to the 
same eigenvalue, with arbitrary coefficients, is also an 
eigenstate with the same eigenvalue. 

Very special properties characterize the antihermitian 
antilinear operators. Within the Dj h or D™ groups only 
the T and V T operators in odd systems belong to such 
a subset (Table [n]). For each antihermitian antilinear 
D™ operator the space of fermion states can be arranged 
in pairs of orthogonal states |4 r (— )}) [ jig] , such 

that, for example, 



T|*(±)) = ±|*(t)>, 



P T |*(±)) = ±|*(T)>. 



(29) 



(30) 



Therefore, the T and V T operators cannot be diagonal- 
ized. In particular, there is no odd fermion state which 
would be invariant with respect to T or P T . 



III. SYMMETRIES OF LOCAL DENSITIES 

Suppose that the Fock-space operator UorW, belong- 
ing to D^ h or D™, respectively, represents a symmetry 
conserved by a mean-field many-particle state | V I'+) or 
|^-), in even or odd fermion spaces, i.e., 



u|*+) = «!*+), 

= u|*_). 



(31a) 
(31b) 



As discussed in Sec. [IE, eigenvalue u can be equal to ±1 
or ±i, and moreover, in odd fermion systems U cannot 
be equal to either T or "P T , i.e, neither the time reversal 
nor the product of inversion and time reversal can be a 
conserved symmetry in odd syst ems. According to con- 
ventions introduced in Sec. II C, in odd systems the hat 
alwa ys de notes one of the operators introduced in 
Sec. ( [I A ), and not one of their partners [Eqs. ( [H|) and 
(p2|)]. Of course, if U is a symmetry of ) then U is 
a symmetry as well, so from the point of view of con- 
served symmetries, any extra study of partner operators 
is unnecessary. 



G 



TABLE III. Symmetry properties of space-spin one-particle operators O belonging to different 
one- dimensional ircoreps of the Dj h and D™ groups. The first column gives names of different 
ircoreps, the second column lists examples of operators O, and the remaining columns give signs in 
the expression U*OU=±d, for operators U given in the column headers. Note that these results 
do not depend on whether operators U act in even or odd spaces, and therefore the Fock-space 
notation is used for them. 



Ircorep 


Space-spin one-particle operators O 


f P 


Rx Ry 


invariants 


x 2 , y A , z z , V£, VjJ, S7i;xiS7 y 6 z ,yiS7 z a x , ziV x a y 


+ + 


+ + 


pseudoinvariants 


xyz, xiV y iV z ,yiV z iVx,ziVxiV y ; iV x a x ,iV y a y ,iV z a z 


+ - 


+ + 


antiinvariants 


xi V x , yiVy , z£V z ; xya z , yza x ,zxa y 


- + 


+ + 


antipseudoinvariants 


xyiV z ,yziV x ,zxiV y ;xa x , ya y , za z 




+ + 


i-covariants 


x; iV ' y a z ,iV z <J y 


+ - 


+ ~ 


i/-covariants 


y\ iV x a z ,iV z a x 


+ - 


- + 


z-covariants 


z; i\7 x a y , iX7 y a x 


+ - 




a;-pseudocovariants 


yz; xiX7 y a y ,iV x ya y , B x yiV y 


+ + 


+ ~ 


t/-pseudocovariants 


xz; yi\7 z a z ,iV y za z ,a y ziX7 z 


+ + 


- + 


z-pseudocovariants 


xy; ziV x a x ,i\7 z xa x , a z xi\7 x 


+ + 




x- ant icovar iant s 


iV x ;ya z , za y 




+ - 


y-anticovariants 


iV y ;xcr z , za x 




- + 


z-anticovariants 


iV z ;xa y , ya x 






x-antipseudocovariants 


yiV z ,ziV y ;a x 


- + 


+ - 


y-antipseudocovariants 


xiV z , zi\7 x ; o y 


- + 


- + 


z-antipseudocovariants 


xiV y ,yiS7 x ; a z 


- + 





Mean-field state I'S) can be characterized by the single- 
particle density matrix p (see Ref. M for the definition), 
for which the symmetry properties (|3l|) imply 



frpU = 



(32) 



independently of eigenvalue u. (Symmetry property ( |32| ) 
does not depend on whether the mean-field state belongs 
to the even or odd fcrmion space, and therefore, we use 
the Fock-space notati on U for the symmetry o perat ors, 
see definitions in Sec. [I A and discussion in Sec. II C.) It 
then follows that the single-particle self-consistent Hamil- 
tonian h[p] is also symmetric with respect to operator U 
0, namely, 



Wh[p]U = h[p}. 



(33) 



Eq. (|3^) implies that if tp is a normalized single-particle 
eigenfunction of h[p], then Up is also a normalized eigen- 
function, both belonging to the same eigenvalue. As a 
consequence, it can be shown [ fl9| that the symmetry is 
preserved during the standard self-consistent iteration, 
provided the entire multiplets of states belonging to the 
same eigenvalue of U are either fully occupied, or fully 
empty. In such a case Eqs. ( p2| ) and fl33|) are fulfilled re- 
peatedly in the successive steps of iteration, and U is a 
self-consistent symmetry. 

Since the one-body density is a fermion-number con- 
serving one-body operator, it can be classified accord- 
ing to one-dimensional ircoreps of D J h or D™ , and this 
can be done both in even and odd systems. This means 



that either the given operator U is a conserved symmetry, 
Eqs. ( pl| ) and (|3^), and the density matrix belongs to the 
given one-dimensional ircorep of the subgroup generated 
by U, or U is a broken symmetry and the density matrix 
has two non-zero components in two different such one- 
dimensional ircoreps. It follows that in odd systems the 
density matrix has always non-zero components in two 
ircoreps corresponding to the time reversal. 

This classification procedure is used below to enumer- 
ate properties of the density matrix when one or more 
Dj h or D™ operators are conserved symmetries. Note 
also, that unlike for the many-body states one does 
not have a freedom to change the phase of the density 
matrix, because it is a hermitian operator independent 
of the phase of the mean-field state it corresponds to. 
Therefore, if the density matrix has non-zero components 
in two ircoreps corresponding to two different eigenvalues 
of an antilinear Dj h or D™ operator, it cannot be trans- 
formed to the form in which it would have been either 
even or odd with respect to this operator. 

A definite symmetry of the density matrix, Eq. (|3^), 
implies certain symmetries for local densities and their 
derivatives. These symmetries are discussed and enu- 
merated in the present section. 

The spin structure of the density matrix is given by 



p(ra,r'cr') = \p(r,r')5 aal 

+ \ s k {r,r') < cr\& k \a' >, 



(34) 



k—x,y,z 



where r=(x,y, z) and r'=(x',y',z') represent three-di- 
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mensional position vectors. When the rotational symme- 
try is preserved one often refers to p(r, r') and Sfc(r, r') as 
the scalar and vector densities, respectively. In our case, 
the rotational symmetry is broken, and we will avoid us- 
ing these terms. Instead, we classify the densities ac- 
cording to the ircoreps of the Dj h or D™ group. As dis- 
cussed above, for the one-body operators only the one- 
dimensional ircoreps are relevant for the classification. 
There are 16 characteristic transformation prope rties of 
the bases for one-dimensional ircoreps. In Table III we 
list all these ircoreps, illustrated by examples of space- 
spin operators of interest, e.g., powers of coordinates, x, 
y, z and gradients, 

Vy, V^. 

The Table also lists explicitly the transformation prop- 
erties of operators belonging to every type of symme- 
try. For example, the sign " — " which appears in row 
denoted by " y-covariants" and column denoted by R x 
means that W x yR x =—y. It can be easily checked that 
the Pauli matrices, a x , & y , a z transform under the sig- 
natures as the x,y,z coordinates, respectively, do not 
change under the inversion, and change their signs under 
the time reversal. Therefore, these can be classified as k- 
antipseudocovariants for k — x,y, z, respectively. Spin- 
dependent operators belonging to other ircoreps can also 
be constructed from the Pauli matrices. Therefore, ex- 
amples of spin-d epe ndent operators are also listed in the 
Table. In Table III we have introduced the same names 



for opera tors as for the bases of one-dimensional ircoreps 
(see Sec. |lIb| ). 

Similarly as in Ref. |^0| , we consider the following local 
densities: 

- particle and spin densities 

p(r)=p(r,r), (35a) 
sfc(r) = s k (r,r), (35b) 

- kinetic and spin-kinetic densities 

r ti W = [V t V' !P (r,r')U„ (36a) 
T Mm {r) = [V fe V' iSm (r, r')] r=r , , (36b) 

- current and spin-current densities 

3k(r) = 1 [V fc - V k )p(r, r')] r=r , , (37a) 
■hi (r) = ^ [(V fc - V k ) Sl (r , r')] r=r , , (37b) 

where each index k, !, or m may refer to either of x, y, or 
z. It follows from the hermiticity of the density matrix p 
that all the above local densities are real functions of r. 
Usually only the traces of kinetic densities, 



T ( r ) = 5Z Tfefc ( r )' 

fc 

T m (r) = ^T fcfcm (r) 



are used in applications. 

When operator U represents a conserved symmetry of 
the density matrix, Eq. (|32|), the transformation rules 
for gradients and spin operators, given in Table 



III, im- 



ply definite transformation rules for the local densities. 
These are listed in Table IV, for all the one-dimensional 



ircoreps of or D™ as indicated in the first col- 
umn. In the second column we show the local densities in 
forms defined by Eqs. (|35|)-(|37|), while the third column 
gives, when possible, the local densities in the traditional 
vector-tensor notation, e.g., 



T = (T x ,T y , T z ), 



J = Jkk, 



(J)ki = -^{Ju + Jik) - -JSki, 



{J)k — ZklmJh 



(39a) 
(39b) 



(40a) 

(40b) 
(40c) 



I in 



Derivatives of densities up to the second order are also 
included in the Table. 

From Table ^ one can read off the symmetry proper- 
ties of various densities. Suppose d(x,y,z) is a generic 
name of one of the densities listed in the second or third 
column, and U is a generic name of one of the Dj h or 
D™ operators listed in the first row. We use the con- 
vention that index % may take any value among x, y or z, 
while indices k ^= I ^= m are arbitrary permutations of x, 
y, and z. If U represents a conserved symmetry, one has 
the following symmetry rule for the density d(x, y, z): 



d{e x x, e y y, e z z) = ed(x, y, z), 



(41) 



where e is the sign listed in Table IV in the row denoted 
by d and column denoted by U . Signs (e x , e y , e 2 ) are 



given in the last row of Table [II, and pertain to two 
Dj h or D™ operators (viz. U and U T ) in two adjacent 
columns. These latter signs give changes of coordinates 
(x, y, z) under the action of U. As the time reversal does 
not affect spatial coordinates, these signs are the same for 
any pair of operators U and U T . One generic Table of 
signs determines, therefore, symmetry properties of any 
local density for any of the Dj h or D™ symmetries being 
preserved. 

For example, symmetry properties of density J xy can 
be found by using indices l—x and m=y (which requires 
k=z) in the row pertaining to fc-covariants. For the con- 
served R z =Rk symmetry we then find in the correspond- 
ing column e=+ and e x =ei=-, e y =e m =—, and e z =e k =+, 



(38a) which gives J xy (-x, -y, z)=J xy (x, y, z). 
(38b) 
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TABLE IV. Symmetry properties of various local densities belonging to different one-dimensional ircoreps of the D^ h or 
D™ groups. For instructions on using the Table see the text explaining Eq. (|4l|). 



T 

Ircorep 


Local densities 


i 


pT 
r 


P, 


H k 


c\ 
Ok 


AT 
"->fc 


til 






AT 




at 


invariants 


P, Tii, ViP, VkJlm 


p,r, Ap, V • J 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ + + 


+ 


pbCLiouiiivai icuito 


T 

"1% 


T 

J 






1 


1 
\ 






1 
\ 


1 






1 1 




antiinvariants 


Tklm, Vjj'j, VfcV;S m 




+ 


— 


+ 




+ 


— 


+ 




+ 


— 


+ - + 


— 


antipseudoinvariants 


ViSi 


V ■ s 




+ 


+ 






+ 


+ 






+ 


+ - - 


+ 


fc-covariants 


Jim, VfeP 


(J)fc, ( J)lm, (Vp)fe 






+ 


+ 










+ 


+ 


- - + 


+ 


fc-pseudocovariants 


Tim, V; V m p, Vfc Jii, 
^ ' iJki, ^ ' iJik 


(VJ) fc , 

(V- J) fc ,(V x J) fc 


+ 


+ 


+ 


+ 


+ 


+ 














k- ant icovar iant s 


jk, V;S m 


(j)fc,(Vxs) fe 




+ 


+ 






+ 




+ 


+ 




- + + 




fc-antipseudocovariants 


Sk , Tiik , Tkii , V ijm , 

V fc V,Si, V?s fc 


(«)*,(T)*,(Vxi)* 
(V(V- S )) fc ,(A s ) fe 


+ 




+ 




+ 






+ 




+ 


- + - 


+ 






(efc, e;, e m ) = 


(- 


--) 


(+- 


--) 


(- 


+ +) 


(- 


+ -) 


(+ 


-+) 


(--+) (+ + 





It is worth noting that symmetry properties (^Tj) which 
correspond to various Dj h or D™ operators, are related 
to one another only by the corresponding group multipli- 
cation rules. Therefore, a specific choice of the conserved 
generators, either for the complete Dj h or D™ groups or 
for any of their subgroups |L4j, leads to a specific set of 
symmetry properties of local densities. 

Symmetry properties (|4l]) can be used for the purpose 
of a continuation of densities from one semi-space into 
the second semi-space, i.e., one can use only space points 
for, e.g., x>0. For two symmetry properties ([fl"|), coming 
from two different symmetry operators (but not from the 
pair U and U T ), one can restrict the space to a quarter- 
space, where two coordinates have definite signs, e.g., 
a;>0 and y>0. Finally, three conserved symmetries allow 
for a restriction to one eighth of the full space with all 
the coordinates having definite signs, e.g., x>0, y>0, and 
y>0. The time-reversal symmetry does not lead to re- 
strictions on the space properties of densities, but, when 
conserved, gives the vanishing of all the antiinvariant, 
antipscudoinvariant, anticovariant and antipseudocovari- 
ant densities, viz., Sk, jk, Tkim for arbitrary fc, I, m as well 
as their derivatives (see Table IV). The possibilities of 



simultaneously conserving one, two, three, or four sym- 
metry operators from the Dj h or D™ groups will be 
discussed in Ref. Q . 

Since density matrix p and single-particle Hamiltonian 
h[p] are always simultaneously invariant under any con- 
served symmetry U, Eqs. ( |32| ) and (|33|), the discussion 
above can be repeated for self-consistent local fields ap- 
pearing in a local mean-field Hamiltonian. Explicit for- 
mulas for symmetry properties of local fields are identi- 



cal to those listed in Table IV, and will not be repeated 



here. In applications, these symmetries appear automat- 
ically when the self-consistent mean fields are calculated 
in terms of densities, cf. Ref |2Ql. 



IV. SYMMETRIES OF SHAPES, CURRENTS, 
AND AVERAGE ANGULAR MOMENTA 



In this section we discuss properties of average values 
of various operators, calculated for the HF many-particle 
state | 'J'). In particular, we consider the electromagnetic 
multipole operators and the total angular momentum; 
the quantities which are used to characterize properties 
of investigated systems. First of all, we enumerate trans- 
formation properties of these operators under the Dj h 
or D j h D operators 



(Sec. Ill 



Similarly as for the density matrix 
the one-body operators discussed in this sec- 
tion belong to one-dimensional ircoreps of Dj h or D™, 
and therefore, their properties do not depend on whether 
the system contains even or odd number of fermions. 



A. Transformation properties of angular momentum 
and multipole operators 

The /c-component of total angular momentum, Ik, 
transforms obviously as fc-antipseudocovariant under Dj h 
or D™, and its transformation rules can be easily read 
off from Table fu. 







TABLE V. Symmetry properties of electric multipole operators Qx^ with respect to operators of 
the Dj h or D™ groups. The results of the symmetry operator WQx^U axe given for three spatial 
directions k—x, y, z. Where applicable, the upper part of the Table gives expressions in terms of 
changed signs of magnetic components, and the lower part gives the equivalent expressions in terms 
of the complex conjugation. 



k 


Rk 




Sk 


si 




X 


(-1) a Qa,- m 


(-l) A+fl QA M 


Qa,- m 


(-1)" 




y 


(-1) A -*Q A ,_ M 


(-i) a Qa m 


(-1)-"Qa,-m 






z 


(-1) m Qa m 


Qa,- m 


(-i) a+m 0a m 


(-1) A 




X 


(-1) A+m Qa m 


(-1) a+m Qa m 


(-1) m 0a m 


(-1)" 




y 


(-i) a Qa„ 


(-i) a Qa„ 


Qa m 


( 


Qx^, 


z 




(-1) m 0a m 


(-1) a+m Qa m 




Qx^ 



For A even (odd), the electric multipole operators 
are even (odd), respectively, under the action of the in- 
version, and are all even with respect to the time reversal, 
i.e., 



(42a) 
(42b) 



The magnetic multipole operators M Xfl have opposite 
transformation properties, 



P^M Xfl P=-{-l) x M Xli , 



(43a) 
(43b) 



Table gives transformation properties pl[ | of Q XfJl 
with respect to operators of the or D™ groups, other 
than T and P. One may note that the electric multipole 
operators are invariant with respect to the Sy symmetry. 
This is of course a consequence of the standard phase 
convention for the rotational irreducible tensor operators 

"HIM 



Qxm = (-i)^Qa,- 

(-1)-"M Ai _, 



M 



IXfj. 



(44a) 
(44b) 



which ensures that the antilinear operator Sy acts as an 
identity upon any irreducible spherical tensor operator. 



B. Average values 

The electric and magnetic moments are defined as 
Qa m = (*|Qa„|*> = / qxn(r) d 3 r, (45a) 
M Xfl = (*|jfir A/ i|*> = [ m Xf ,(r)d 3 r, (45b) 



where "J) is a many-body mean-field state, and q Xtl (r) 
and m x ^(r) are the corresponding moment densities: 



q\»(r) = ep(r)Q XjJi (r), 
m Xfl (r) = n N 2J (5sSfcV fc QA M (r) 



(46a) 



k=x,y,z 



j^l9dk{r x VQ A(I (r)), 



(46b) 



and e, and <7; are the elementary charge, and the spin 
and orbital gyromagnetic factors, respectively In def- 
initions (f46[), multipole functions plj (solid harmonics) 
have the standard form: Q Xtl (r)—P : Y Xtl (0, </>). 

Similarly, the mean value of the fc-component of total 
angular momentum (in units of h) reads 

h = <*|4|*> = / (e klm r lJm (r) + ±s fe (r)) d 3 r (47) 



We may now co mbi ne symmetry properties of densities 
p, s, and j, Table III, with those of multipole operators, 
Table [v| to obtain symmetry conditions obeyed by the 
electric and magnetic moments, and by the average an- 
gular momenta, for given conserved symmetries of the 
D^h or D™ groups. In doing so, we have to remember 
that since the electric multipole operators are time-even, 
the corresponding electric moments depend only on the 
time -even component of the density matrix, as given in 
Eq. (46a). This is so irrespective of whether the time re- 
versal is, or is not a conserved symmetry, or whether the 
system contains even or odd number of fermions. There- 
fore, the time reversal does not impose any condition on 
the electric multipole moments. On the other hand, with 
the time-reversal symmetry conserved, which may occur 
only for even systems, all magnetic moments and aver- 
age angular momenta must vanish, because they depend 
only on th e time-odd component of the density matrix, 
Eqs. (ggbh and @. 
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TABLE VI. Conditions fulfilled by the electric and magnetic multipole moments, Qx^ and 
M\^, for conserved Dj h or D™ operators. Where applicable, the upper part of the Table gives 
expressions in terms of changed signs of magnetic components, and the lower part gives equivalent 
expressions in terms of the complex conjugation. 



h. 

k 














c 






AT 

6,,. 


X 


Q\n— 


(-i) a <3a,-p 


Qa m — 


(-1) A Qa,- 


M 


Qa m — 




M 




Qa,-m 


X 


fl/f 


(,-lJ Mx,-f, 




/ 1 \A 71^ 

— (, — lj MX,- 


-M 


7\yf 


71 4" 


_ M 


71 /f 


A^A,-m 


y 


Qam=(- 




Qxij,~ 


(-1) a -"Qa,- 


M 


Qa m = 


(-1)-"Qa,- 


M 


Qxij.= 


(-1)-"Qa,- m 


y 




1) M M Ai _^ 


71 ^ 

-Ma m = 


/ 1 \A — U 71/f 

-(-1) A4a,- 




71 ^ 

JWa m = - 


/ 1 \ — A 4 71/f 

-(-1) ^A1a,- 




71 //" 

A^A^ = 


(-1) M M A ,_ M 


z 


Q\n= 


(-1)"0a„ 


QxiJ.— 


(-1) m Qa„ 




Qa m = 


(-1) A+m Qa m 




Qa m =( 


-1) A+m Qa m 


z 




(-1) m M Am 


Ma m = 


-(-1)"Ma, 




Af AM =- 


(_1)^+m Mam 




M Am =( 


-1) A+m Ma m 


X 


Qa m = (- 


1) A+m Qa„ 


Qa m = 


(-1) a+ "Qa„ 




Qa m = 


(-1) m Qa m 




Qa m = 


(-1) m 0a„ 


X 


M Am =(- 


1) A+M M A * M 


Ma m = 


-(-l) A+fl M A * M 




Af AM = 


-(-1)^ 




M Am = 


(-l) M AfA* M 


y 


Q\n= 


C 1) Q ambda/j.. 


Qa m = 


(-1) A Qa m 




Qa m = 


Qa/j 




Qa m = 


Qa/j 


y 






Ma m = 


-(-1) A A^ M 




Af AM = 






M Am = 


Ma* m 


z 


Q\n= 


(-1) m Qa m 


QxiJ.~ 


(-1) m Qa m 




<3a m = 


M) A+m Qa m 




Qa m =( 


-1) A+m Qa m 


z 


M Am = 


(-1)^Ma m 


M Am = 


-(-1)"Ma m 






(_1)A+/- Mam 




M Am =( 


-1) a +^Ma m 



For the conserved parity, one obtains the standard con- 
ditions: 



Qxtj, 



-(-1) a M Am , 



(48a) 
(48b) 



i.e., odd electric and even magnetic moments must van- 
ish. Similar symmetry properties with respect to other 
symmetries of the Dj h or D™ groups are collected in 
Table 0. 

Within the standard phase convention of Eq. (|44|) , only 
a conservation of the y-T-simplex symmetry, Sy , enforces 
the reality of all multipole electric and magnetic mo- 
ments. In such a case, the lower part of Table VI gives 



at a glance all the multipole moments which must vanish 
whenever any other symmetry is additionally conserved; 
these are those for which the phase factors are negative. 
On the other hand, a conservation of the x-T-simplex 
symmetry, S*J , enforces the equality of negative and pos- 
itive magnetic components. In this case, a conservation 
of any additional symmetry puts to zero the multipole 
moments with negative phase factors appearing in the 
upper part of the Table. Of course, numerous other com- 
binations of conserved symmetries can be considered, for 
example, a conservation of the y-simplex symmetry, S y , 
gives real electric moments and imaginary magnetic mo- 
ments. 



Since conditions listed in Table VI depend only on the 
parity of A and on the parity of /x, and since condition 
(||) allows us to consider only non-negative values of 
/i, one has only six types of the symmetry properties 
of multipole moments with respect to the Dj h or D™ 



operators. These six types are listed in Table VII for 
electric and magnetic moments. Column denoted by the 
identity operator E gives the properties resulting solely 
from condition (44), while the remaining columns give 
properties of moments when one of the non-identity Dj h 
or D™ operators is conserved. 

In the same Table we also give symmetry properties 
of the Cartesian components of the average angular mo- 
menta Ik (p7|). Although the symmetry properties of the 
angular momentum are identical to those of the dipole 
magnetic moment, explicit values shown for its Cartesian 
components allow for a simple visualization of a direc- 
tion taken by the angular-momentum vector when vari- 
ous Dj h or D™ operators are conserved. In particular, 
one can see that a conservation of any of the signature or 
simplex operators for a given axis enforces the angular- 
momentum direction along that axis, while a conserva- 
tion of any T-signature or T-simplex operators allows for 
a tilted angular momentum in the plane perpendicular 
to the given axis, cf. Ref. |23| ]. On the other hand, none 
of these operators may be conserved if the angular mo- 
mentum is to be tilted beyond any of the x-y, y-z, and 
z-x planes. Note, however, that the above tilting con- 
ditions pertain to the reference frame, and not to the 
principal axes of the mass distribution. An appropriate 
choice of the reference frame, as discussed below, has to 
be performed in order to relate the conserved or D™ 
operators to the direction of Ik with respect to the mass 
principal axes. 
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TABLE VII. Properties of electric multipole moments Qx^, magnetic multipole moments Mx^, 
and average angular momenta, Ik for conserved Dj h or D™ operators. Symbols C, R, I, or denote 
values which can be, in general, complex, real, imaginary, or zero, respectively. 
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Independently of any Dj h or D™ symmetry break- 
ing, the reference frame in the space coordinates can be 
chosen in such way that some of the moments have sim- 
ple forms. For example, a shift of the reference frame 
can bring all electric dipole moments to zero (this corre- 
sponds to using the center-of-mass system of reference), 



i.e. 



Qw = KQii - 9fQu = 0. 



(49) 



Similarly, a suitable rotation of the reference frame can 
bring the electric quadrupole moments Q2fi to the prin- 
cipal axes, where 



= 3Q 2 i = QQ22 = 0. 



(50) 



On the other hand, for some conserved symmetries, 
these conditions can be automatically satisfied. For ex- 
ample, conservation of the T>2h group (i.e., simultane- 
ous invariance with respect to operators P, R x , and 
R y ) ensures that the center-of-mass (fl9|) and principal- 
axes ([50]) conditions are automatically satisfied, see Table 
VII. Therefore, the breaking of the Dj h or D™ symme- 
tries may have non-trivial physical consequences only for 
higher electric moments; starting from Q30, if the parity 
is broken, or starting from Q41, for example, if the parity 
is conserved. In other words, the Dj h or D™ symmetry 
breaking will not lead to new classes of low-multipolarity 
shapes. Nevertheless, such symmetry breaking will im- 
mediately be reflected in values of magnetic moments, 
whenever the time reversal is broken too. 



V. CONCLUSIONS 

In the present study we have presented applications 
of point groups based on the three mutually perpendic- 
ular symmetry axes of the second order, inversion, and 
time reversal, to nuclear structure problems. We have 
discussed properties of the corresponding single Dj h and 
double groups in describing even and odd fermion 
systems, respectively. We have enumerated their repre- 
sentations, both for many-body states and for the single- 
particle operators, and reviewed properties of group op- 
erators when they are represented in the fermion Fock 
space. 

Consequences of conserving individual Dj h or D™ 
symmetries have been enumerated for: (i) space sym- 
metries of local one-body densities, (ii) electric and mag- 
netic multipole moments, and (iii) average values of the 
angular-momentum operators. This gives information 
about the nuclear shapes and matter-flow currents in 
states obeying one or more of the Dj h or D™ symme- 
tries, and allows for selecting appropriate conserved sym- 
metries in descriptions aiming at various physical phe- 
nomena. 
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APPENDIX: 

In this Appendix we explicitly construct irreducible 
representations of the E)JP group by using the example 
of the HO basis (Sec. |II D| ), and we illustrate the Wigner 
classification of groups that contain antilinear operators 
(see Chap. 26 of Ref. J8|). The results of such an analysis 
were used in Sec. |J. 

We consider here only the two-dimensional spinor rep- 
resentations, appropriate for the odd-fermion systems 
and in particular for the single-particle states. From 
Eqs. (p7) one finds representation matrices T(JJ) (where 



operators £YeDJ? h of Sec. II C form the double group D^ h ), 
in the two-dimensional invariant subspace spanned by 



i) and \n x n y n z , . 



\). We have 



-T(S) = a , 



T(V 

T(K k ) = -T(TZ k ) = -i{-\) Nk a k 
Y{S k ) = -T(S k ) = -i(-l) n *a k , 



Co, 



(Ala) 
(Alb) 
(Ale) 
(Aid) 



where <7o is the identity 2x2 matrix, a k for k — x,y, z are 
the standard Pauli matrices, and symbols N x , N y , and 
N z refer to n y +n z , n x +n z , and n x +n y , respectively. 
The characters of the classes are 

X {£) = -X{£) = 2, (A2a) 
x {p) = - x (p) = 2{-l) n * +n » +n *, (A2b) 
x({6»,£*}) = X({S X ,S X }) = 0, (A2c) 

x({n y ,n y }) = x({Sy,s y }) = o, (A2d) 

X({n z ,n z }) = x({Sz,S,}) = O. (A2e) 

One can see, that only the characters of V and V depend 
on quantum numbers n x , n y , and n z that define the in- 
variant subspaces; more precisely, they depend only on 
the parity of the sum n x + n y + n z , i.e., on the total parity 
of basis states. Therefore, the only two spinor represen- 
tations of D™ can be labeled by the eigenvalues of the 
parity operator V . Let us also note that all characters 
are real. 

If we introduce the time reversal, T, into the ensem- 
ble of the linear operators belonging to D^ h we obtain 
the D™ group with 16 new antilinear elements IA T =UT , 
Sec. [I C . To study properties of the representations of 
the group, one has to consider representations pro- 
vided by matrices 



f (ti) = Y{A- 1 UA)*, 



(A3) 



where A is one of the antilinear elements of D™ (see 
Q). It is most convenient to take T itself as A; we then 
have simply 



t{u) = r(uy 



(A4) 



as T commutes with all WeD^ h . In such a case, matrices 
are just complex conjugates of T(W), and therefore 
the characters of representation T are exactly the same 
as those of T, because they are all real, see Eqs. (|A^). 
Therefore these two representations are equivalent, and 
a matrix (3 exists which brings by a similarity transfor- 
mation all matrices T(U) to T{U)=T(U)*, 



U e D°. 



(A5) 



Now, as shown by Wigner |£j], there are only two cases 
possible: either 



1313* = +r(f 2 



(3p* 



-r(T 2 



(A6) 



(A7) 



Matrix (3 can easily be found fro m the explicit expressions 
for matrices T(U) given in Eqs. (Al), and it reads 



-1 

1 o 



—ze r a 7 , 



(A8) 



Choosing the phase factor e^^l in ( [A§| ) is equivalent to 
a change of phase of the \n x n y n z , s z = + ^ ) sta tes, and to 
a change in the phase convention in Eq. (27b). 

It is easy to demonstrate that with this form of the 
matrix (3, Eq. (|A6|), and not (A7) holds. In Wigner's 
classification this case leads to what is called the corep- 
resentations of the "first kind" : any representation T of 
the group D^ h can be completed to a corepresentation of 
the full D™ group by defining 



T{U T ) = T{U)(3. 



(A9) 



Note that taking U=£ one gets T(T)—(3, so (3 is, of 
course, just the matrix representing T itself. 

After Wigner, the term corepresentation is used here 
because the representations of groups containing antilin- 
ear operators are not representations in the usual sense. 
To see this, let us consider an orthonormal set of states 
{|<?!>i)} constituting a basis of a representation T. Let U 
be any linear, and U 1 any element of the group. Then, 
because of antilinearity of U T one has 

{u T -u')^) = ]Tw T r(w%|0 4 ) = ^Tr(w%-zF|0 4 ) 

i i 

= Y,nu%nu T )^ k ) 



ik 



[r(u T )r(u'r 



kj 



and, consequently, 

T{U T -U') =T(U T )T(U'y, 
to be compared with 



(A10) 



(All) 
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r(w -u') = Y{ii)Y{u'), 



(A12) 



which holds for the "usual" representations. The pres- 
ence of c omplex conjugation on the right-hand side of 
Eq. (All) implies that the homomorphism between the 
group multiplication and the multiplication of represen- 
tation matrices no longer holds when the group contains 
antilinear operators. This is not surprising in view of the 
fact that matrices, by construction, always act on vectors 
(columns of numbers) linearly. 

In conclusion, there are only two spinor corepresenta- 
tions of D™, and they can be labeled, as is also the case 
for the D?? h group, by one quantum number only (parity). 



[19] 
[20] 

[21] 



[22] 



[23] 
[24] 



to the one used by Bohr and Mottelson H|. 
G. Ripka, Adv. Nucl. Phys. 1, 183 (1968). 
Y.M. Engel, D.M. Brink, K. Goeke, S.J. Krieger, and D. 
Vautherin, Nucl. Phys. A249, 215 (1975). 
D.A Varshalovitch, A.N. Moskalev, and V.K. Kersonskii, 
Quantum Theory of Angular Momentum (World Scien- 
tific, Singapore, 1988). 

We adopt the general phase convention ]2lj valid for arbi- 
trary irreducible tensor operators. However, note that for 
the integer-angular-momentum irreps, the signs of mag- 
netic components appearing in the phase factors can be 
arbitrarily changed. 

S. Frauendorf, to appear in Reviews of Modern Physics. 
A. Bohr and B.R. Mottelson Nuclear Structure (Vol. I) 
(W.A. Benjamin, Inc., 1969), p. 19. 



P. Ring and P. Schuck, The Nuclear Many-Body Problem 
( Springer- Verlag, Berlin, 1980). 

M. Hamermesh, Group Theory (Addison- Wesley, Read- 
ing, Mass., 1962). 

G.F. Koster, J.O. Dimmock, R.G. Wheeler and H. Statz, 

Properties of the Thirty-Two Point Groups (M.I.T. 

Press, Cambridge, Massachusetts, 1963). 

L.D. Landau and E.M. Lifschitz, Quantum Mechanics 

(Non-relativistic Theory) (Pergamon Press, 1981). 

J.F. Cornwell, Group Theory in Physics (Academic 

Press, 1984). 

Ch.W. Curtis, I. Rainer, Representation Theory of Finite 
Groups and Associative Algebras (Interscience Publish- 
ers, New York, 1962). 

C.J. Bradley, A.P. Cracknell, The Mathematical Theory 
of Symmetry in Solids (Clarendon Press, Oxford, 1972). 
E.P. Wigner, Group Theory and Its Application to the 
Quantum Mechanics of Atomic Spectra (Academic Press, 
New York, 1959). 

P. Bonche, H. Flocard, and P.-H. Heenen, Nucl. Phys. 
A467, 115-35 (1990). 

J. Dobaczewski and J. Dudek, Comp. Phys. Commun. 
102, 166 (1997); 102, 183 (1997). 

J.W. Negele, Lecture Notes in Physics 40 (Springer, 
Berlin, 1975) pp. 285, 288. 

R.M. Dreizler and E.K.U. Gross, Density Functional 
Theory (Springer, Berlin, 1990). 

P. Quentin and H. Flocard, Annu. Rev. Nucl. Part. Sci. 
28, 523 (1978). 

J. Dobaczewski, J. Dudek, S.G. Rohozihski, and T.R. 
Werner, the following paper. 

A. Messiah, Quantum Mechanics (Wiley, New York, 
1962). 

In this article, superscript T always means that the cor- 
responding operator contains factor T; this should not 
be confused with the notation for the transposition op- 
eration. 

M. Girod and B. Grammaticos, Phys. Rev. C27, 2317 
(1983). 



The phase convention implied by Eq. (27b) agrees with 
that of Varshalowitch et al. |H| , and thus it is opposite 



14 



